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The FLIC Overlap Quark Propagator 
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FLIC overlap fermions are a variant of the standard (Wilson) overlap action, with the FLIC (Fat 
Link Irrelevant Clover) action as the overlap kernel rather than the Wilson action. The structure 
of the FLIC overlap fermion propagator in momentum space is studied, and a comparison against 
previous studies of the Wilson overlap propagator in quenched QCD is performed. To explore the 
scaling properties of the propagator for the two actions, numerical calculations are performed in 
Landau Gauge across three lattices with different lattice spacing a and similar physical volumes. 
We find that at light quark masses the actions agree in both the infrared and the ultraviolet, but at 
heavier masses some disagreement in the ultraviolet appears. This is attributed to the two actions 
having different discretisation errors with the FLIC overlap providing superior performance in this 
regime. Both actions scale reasonably, but some scaling violations are observed. 

PACS numbers: 12.38.Gc, 11.15.Ha, 12.38.Aw, 14.65.-q 



I. INTRODUCTION 

The quark propagator in momentum space is a funda- 
mental quantity in QCD. Its infrared structure provides 
insight into the dynamical generation of mass due to the 
spontaneous breaking of chiral symmetry within QCD. 
Its ultraviolet behaviour allows one to obtain the running 
quark mass 0, 0, and OPE condensates . Lattice 
QCD allows a direct probe of the non-perturbative quark 
propagator. The momentum-space quark propagator has 
been studied previously usingdifferent gauge fixing and 
fermion actions [3J, |5|, |fj, |7], |8|, |9j, [lfj , including studies of 
the (standard) Wilson overlap in Landau gauge |3.I^.ITo|. 
The FLIC overlap [HI has been shown to possess com- 
putational advantages over the Wilson overlap, and we 
are interested in comparing the (quenched) Wilson over- 
lap and FLIC overlap propagators to see if there is any 
difference in their scaling properties. 

In the continuum, the tree-level (A tl (x) — 0) quark 
propagator is identified with the (Euclidean space) 
fermionic Greens function, 



< + m°)A f (x,y) = 5 4 {x-y), 



(1) 



where m° is the bare quark mass. In momentum space 
this equation is solved straightforwardly, 
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Denote S^(p) = Af(p) to be the tree-level propagator 
in momentum space. Then in the presence of gauge field 
interactions, define Share (p) to be the (fourier transform) 
of the (interacting) fermionic Green's function, 

(p + m°)A h l ™ c {x,y) = 5\x-y). (3) 

We define the mass function M(p) and the bare renor- 
malisation function Z{p) such that the bare quark prop- 
agator has the form 
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Then for the renormalisation point the renormalised 
quark propagator is given by 



S((P) = 



ij> + M(p) 



= Z 2 (C,a)S baxe (p), 



(5) 



where Zq(p) is the (^-dependent) renormalisation func- 
tion, and Z2(Ci a ) is the wave function renormalisation 
constant, which depends on £ and the regulator param- 
eter a. The a-dependence of Share i s implicit. Z2(C,ci) is 
chosen such that 



z c (0 = i. 



(6) 



As S^(p) is multiplicatively renormalisable, all of the in- 
dependence is contained within Z^(p), that is, the mass 
function is £- independent. 



II. OVERLAP FERMIONS 

The overlap formalism [l2l Il3l ITS . fl5| in the vector- 
like case leads to the following definition of the massless 
overlap-Dirac operator [l6j. 

1 -(l+ 75 e(i/)). (7) 
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Here, e(H) is the matrix sign function applied to the over- 
lap kernel H. The kernel can be any Hermitian version 
of the Dirac operator which represents a single fermion 
species of large negative mass. The standard choice is 
the Hermitian Wilson-Dirac operator (setting a = 1), 



H w =7s(y + ^A-m w ), 



(8) 



where y is the central covariant finite difference operator, 
and A is the lattice Laplacian, or Wilson term. In terms 
of the parallel forward (backward) transport operators 

T^(x) = Ua{x)i>{x + e„), (9) 

tfil>(x) = Ut(x-e lt )il>(x-e ll ), (10) 
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we have 



III. THE OVERLAP PROPAGATOR 



2 C^p 



(ii) 

(12) 



It has been shown that alternative kernel choices, in par- 
ticular the FLIC (Fat Link Irrelevant Clover) action , 
can accelerate the evaluation of the overlap-Dirac oper- 
ator fill ITsf . The FLIC fermion action is a variant of 
the clover action where the irrelevant o per ators are con- 
structed using APE-smeared links [H Hj El HI, and 
mean field improvement |23j is performed. The Hermi- 
tian FLIC operator is given by 

tffiic = T5(y mfi + i(A« fi - la • i£ fi ) - m w ), (13) 

where the presence of fat (smeared) links and/or mean 
field improvement has been indicated by the super- and 



subscripts. In this work we choose t7 M „ 
use a standard one-loop F^ U1 



Wl^lA and 
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(14) 



+ U Vj Jx) + U_,^(x)), (15) 



where U^. v (x) is the elementary plaquette in the +/x, +f 
direction. 

The APE smeared links U^(x) constructed from U^x) 
by performing n smearing sweeps, where in each sweep 
we first perform an APE blocking step, 



#(x) = (1-q)_+ 



(16) 



followed by a projection back into SU(3),ujf\x) = 

"Piy^ {x))- In this work, the projection is performed by 
choosing the matrix which maximises the following gauge 
invariant measure, 



Ujf\x) 



max ReTrtU'JxW^Ux)). (17) 
U'esu(3) M ' 



As it is only the product not that matters |2J], we fix 
a = 0.7 and only vary n. 

Mean field improvement is performed by making the 
replacements 



U^x) 



u 



U«(x) 



K{x) 



(18) 



where uq and Uq are the mean links for the standard and 
smeared gauge fields. We calculate the mean link via the 
fourth root of the average plaquette, 



Uo 



(-ReTrU^(x))l^. 



(19) 



It is easily seen that the continuum massless quark 
propagator anti-commutes with 75, 



{75,<Sb are (p)| mO=0 } = °' 



(20) 



A straightforward consequence of the Ginsparg- Wilson 
relation is the inverse of the overlap Dirac operator sat- 
isfies 



{ 75 ,£>- 1 } = 2 75 



(21) 



Using either a Wilson or FLIC kernel, and noting that 
the mean link is a function of a and that uo(a), Uq(o) — » 1 
as a — > 0, we obtain 



lim D — 

a^o 2m, 



■p. 



(22) 



It is then natural to define the (external) massless bare 
overlap propagator on the lattice as 0, l2fij 



Share (p)|m°=0 

as we then have that 



2m, 
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{7S,Share(p)| m o =0 } = °> 



(23) 



(24) 



as in the continuum case. The massive overlap Dirac 
operator is given by |26| 



D (n) = (l-n)D +n, 



(25) 



with \fj,\ < 1 representing fermions of mass oc The 
massive (external) bare overlap propagator is defined 
as _25] 



Share (P) 



1 



2m w (l — n) 
and with the identification 



2m v 



satisfies 



bare 



(p)=5bate(P)Lo =0 +m°. 



(26) 



(27) 



(28) 



In order to construct M (p) and Z(p) on the lattice, we 
first define £>(p),C M (p) such that 



Shared) = -ip(p)+B( P ). 



Then 



CM 
B{p) 



n s n c 



Tr^Sbareb)], 

Tr[S barc (p)], 



(29) 

(30) 
(31) 
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where the trace is over colour and spinor indices only, 
and n s , n c specify the dimension of the spinor and colour 
vector spaces. Now, define the functions B(p), C M (p) such 
that the inverse of the bare (lattice) quark propagator has 
the form 



B(p). 



Then it is easily seen that 

C = ^ 

C 2 + B 2 ' 



B = 



B 



C 2 + B 2 '' 



(32) 



(33) 



where C 2 =C-C. 

The kinematical lattice momentum 
that at tree-level 



is defined such 



(34) 



that is q^p) = C^ 0) (p),to° = B<°)(p). We note that the 
simple form of these relations is one of the advantages 
of the overlap, as the absence of additive mass renor- 
malisation prevents the need for having to perform any 
tree-level correction (necessary most other fermions 
actions 0,0, 0)i outside of identifying the correct mo- 
mentum variable q. Now, we define A(p) such that 



S*L(p)=iiA(p) + B(p). 



(35) 



The mass function M(p) and renormalisation function 
Z(p) may then be straightforwardly constructed, 



Zip) 



AipY 



M(p) 



Bjp) 
AipY 



(36) 



Volume 




a (fm) 




^apc 


^Mow 


Phys. Vol. (fm 4 ) 


8 3 x 16 


4.286 


0.200 


0.8721 


7 


20 


1.6 3 x 3.2 


12 3 x 24 


4.60 


0.120 


0.8888 


4 


15 


1.44 3 x 2.88 


16 3 x 32 


4.80 


0.096 


0.8966 


4 


5 


1.536 3 x 3.072 



TABLE I: Parameters for the three different lattices. All use a 
tadpole improved Luscher Weisz gluon action. Lattice spac- 
ings are set via a string tension analysis incorporating the 
lattice Coulomb term. Shown is the lattice volume, gauge 
coupling /3, lattice spacing, the mean link, the number of APE 
sweeps, projected modes and the physical volume. All lattices 
have m w = 1.4. 



choice of m w = 1.4 for the FLIC overlap was slightly 
different from that chosen for the Wilson overlap). The 
tree level propagator is calculated directly by setting the 
links and the mean link to unity. The kinematical lat- 
tice momentum q is obtained numerically from the tree 
level propagator, although it could equally well have been 
obtained from the analytic form for q derived in Ref. |9(. 

Each lattice ensemble consists of 50 configurations. 
The lattice Green's function is obtained by inverting the 
FLIC overlap Dirac operator on each configuration using 
a multi-mass CG inverter |30| . The Fourier transform is 
taken to convert to momentum space, and the the bare 
quark propagator is obtained from the ensemble aver- 
age. The quark propagator is calculated for 15 different 
masses. The details of the FLIC overlap parameters used 
are presented in Table ILTI 



IV. SIMULATION DETAILS 

While the FLIC overlap and Wilson overlap actions are 
both free of 0(a) errors (at both zero and finite quark 
mass), the actions may in general differ at 0{a 2 ), and 
as such there may be some difference in their properties 
at finite a. We conduct a comparison between the two 
actions by calculating the FLIC overlap propagator on 
the same set of quenched (Luscher- Weisz) lattices as in 
the Wilson overlap studies 0, . All the lattices have 
approximately the same volume, and the details of all 
three lattices are given in Table [I] 

Landau gauge is chosen for the gauge fixing. An im- 
proved gauge fixing scheme [27} is used, and a Conjugate 
Gradient Fourier Acceleration algorithm is chosen 
to perform the gauge fixing. We use periodic boundary 
conditions in the spatial directions and anti-periodic in 
the time direction. The matrix sign function is evalu- 
ated using the Zolotarev rational polynomial approxima- 
tion [23 , of degree (typically 8) chosen to give an accu- 
racy of 2.0 x 10~ 8 within the spectral range of the kernel, 
after projecting out low-lying modes. 

At tree level the FLIC overlap and Wilson over- 
lap propagators have identical behaviour (although the 





P = 4.286 


P = 4.60 


P = 4.80 


m° 


Mi 


0.00622 


0.00400 


0.00305 


18 


P-2 


0.01244 


0.00800 


0.00610 


36 


H3 


0.01866 


0.01200 


0.00915 


54 


/'4 


0.02488 


0.01600 


0.01220 


72 




0.03110 


0.02000 


0.01525 


90 


jUe 


0.03732 


0.02400 


0.01830 


108 


/'7 


0.04354 


0.02800 


0.02134 


126 


/'« 


0.04976 


0.03200 


0.02439 


144 


/■'<) 


0.06220 


0.04000 


0.03049 


181 


M10 


0.07464 


0.04800 


0.03659 


217 




0.09330 


0.06000 


0.04574 


271 


fJ-12 


0.12439 


0.08000 


0.06098 


362 


j«13 


0.15549 


0.10000 


0.07623 


452 


M14 


0.18659 


0.12000 


0.09148 


543 


fJ-15 


0.21769 


0.14000 


0.10672 


633 



TABLE II: For each P, displayed is the 15 lattice masses \x. 
The final column shows the bare quark mass in MeV, which 
is approximately the same for each lattice (see Eq. 1271 '). 
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V. RESULTS 

We now turn to the results of our investigation. The 
mass function M(p) and the renormalisation function 
Z((p) are calculated on each of the lattices. A cylinder 
cut 9] is applied to all the data, to reduce the effects of 
rotational symmetry violation. The cylinder cut selects 
momenta with Ap < where rii is the spatial extent 
of the lattice, and Ap = |p|sin6*p is a measure of the 
distance from the hypercubic diagonal n = =(1, 1,1,1). 
The angle 9 p is determined by \p\ cos 9 = p ■ h. The full 
results for each of the lattices are displayed in Fig. ^ for 
M(p) and Fig.[21for Z^(p). The renormalisation point £ 
for (p) is chosen to be approximately p = 3 GeV and 
q = 5.5 GeV when shown against q. 

While the quark masses obtained from the FLIC over- 
lap and Wilson overlap calculations are similar, they are 
not matched. In particular, the lightest five FLIC overlap 
masses are significantly lighter than the corresponding 
Wilson overlap masses. In order to compare the two ac- 
tions, for each fixed momenta value we perform a quartic 
fit of both M(p) and Z(p) as a function of m° = 2m w /z, 
(excluding the lightest two masses for the FLIC Overlap, 
so that the lightest bare mass included in the fit is ap- 
proximately the same for both actions). The data and 
fits for the smallest 10 momenta on each of the lattice 
are shown in Fig. [31 for M(p) and Fig. 0] for Z(p). 

We see that for the curves which correspond to larger p 
the behaviour is essentially linear in the bare mass. How- 
ever, for the curves which correspond to smaller p we see 
that at heavy quark masses the behaviour is linear, while 
at moderate to light quark masses there is an increasing 
amount of curvature present. This is consistent with the 
behaviour of the hadron spectrum predicted by chiral 
perturbation theory. However, chiral perturbation the- 
ory does not give us a prediction of QCD type quantities 
such as M(p). In this work we used a (quartic) polyno- 
mial fit function, which would not take into account any 
non-analytic behaviour that might appear near the chiral 
limit. Using a reliable formula for the extrapolation of 
M(p) and Z(p) requires further study, as we must disen- 
tangle the effects of finite volume, quenched chiral logs 
and genuine curvature. We note at this point that on 
the P = 4.286 lattice the Wilson overlap masses change 
order at large quark mass, an indication that the FLIC 
overlap is superior for heavy quarks. 

Having obtained the fit functions M(p,m°) and 
Z(p,m°) we then choose m° = 0.0,0.05,0.1,0.2,0.4,0.6 
GeV and compare the functions at these values of the 
bare mass. In this way we can compare the two actions 
with matched input masses. Results are shown in Figure 
[3 plotted against the kinematical momenta q. 

First we examine the mass function, plotted vs p (as 
in [U). We see that the actions agree well in the infrared, 
and are both flat in the ultraviolet. At the lighter quark 
masses the two actions agree across the whole momen- 
tum spectrum, indicating the renormalisation of the bare 
mass is similar for the two actions. This is further sup- 



ported by comparing the infrared behavior of the mass 
functions for the two actions at all quark masses. 

For the heavier quark masses different asymptotic val- 
ues of the quark mass are observed. The difference be- 
comes more pronounced as the quark mass increases. 
What this shows is that with a matched input bare mass, 
the low energy properties of the two lattice actions are 
the same. That is, where the lattice approximation is 
good, the actions agree. There are well known difficul- 
ties associated with putting heavy quarks on the lattice, 
as the Compton wavelength of the quarks becomes small. 
It is in this area that we are likely to see sensitivity to the 
discretisation errors, and this is shown by the data, par- 
ticularly in the shape of the momentum dependence. In 
the very infrared the two actions agree for all masses, but 
at heavy quark mass as the momenta increases they begin 
to differ, with the coarser lattice showing a more rapid 
divergence (as expected). Given the relatively strange 
behaviour of the Wilson overlap mass function on the 
coarsest lattice at the heaviest masses, we can conclude 
that the FLIC overlap inherits the benefits of the im- 
proved ultraviolet behaviour of the FLIC kernel, due to 
the APE smearing of the irrelevant operators. 

The comparison of the renormalisation function (which 
we plot vs q) between the two actions leads to similar con- 
clusions as illustrated in Fig. |S] Two things should be 
noted however. The first is that the two actions disagree 
in the ultraviolet regardless of the quark mass. This sug- 
gests that Z(p) is more sensitive to the cutoff than M(p). 
Secondly, we note that we are comparing Z(p) and not 
Z^(p). There is a good reason for this, as by fixing the 
behaviour of Z^(j>) at high momenta and then compar- 
ing the two actions, we would have been led to believe 
(spuriously) that the actions agreed in the ultraviolet but 
not in the infrared. This also tells us that if we choose 
C in the ultraviolet, then ^(Ci a ) 1S different for the two 
actions. 

Using the same technique, we can test the scaling of 
the mass function and renormalisation function. The 
results for both the FLIC overlap and Wilson overlap 
for the three quenched lattices, calculated at mP = 
0.0,0.1,0.2,0.4,0.6 GeV, are displayed in Fig. © Ex- 
amining the FLIC Overlap mass function first, we see 
that the scaling is best at light to medium quark mass. 
At heavy masses the two finer lattices agree reasonably 
well, but the coarsest lattice shows scaling violations. At 
the lighter quark masses the agreement between the lat- 
tices is best in the ultraviolet and the very infrared. In 
the intermediate momentum regime some difference in 
M(p) between the lattices appears. The difference seems 
to increase as the quark mass decreases. One possible 
cause for this difference is in finite volume effects, as the 
j3 = 4.60 lattice had a smaller spatial volume than the 
(3 = 4.80 lattice. If the smaller volume resulted in an en- 
hancement of the mass function, then this would explain 
the difference between the finer lattices at lighter masses, 
as the lighter quark masses would be more sensitive to 
the small volume. 
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The Wilson overlap mass function displays similar scal- 
ing properties to the FLIC overlap mass function, but the 
disagreement between the lattices is in general larger, es- 
pecially at heavy quark mass. We note that at the lighter 
quark masses, the FLIC Overlap shows noticably better 
scaling between 2 < p < 3 Gev, while the Wilson Over- 
lap has slightly better scaling for 1 < p < 2 GeV. The 
Wilson Overlap also appears to scale better in the chiral 
limit. This could be due to differing finite volume effects. 
Another, more likely potential source for the difference is 
in the chiral fits to the data. As mentioned earlier, we use 
a polynomial fit function. This may be a source of sys- 
tematic error at light quark mass, as any non-analyticity 
that might appear in the mass function as we approach 
the chiral limit would be unaccounted for. 

The FLIC overlap renormalisation function Z(q) shows 
agreement between all three lattices in the ultraviolet, 
but in the infrared there are differences between the three 
lattices. Once again, these differences are largest in the 
intermediate momentum regime. In the very infrared the 
disagreement is somewhat smaller, in particular at the 
lighter quark masses. The Wilson overlap renormalisa- 
tion function again shows similar behaviour, but in the 
intermediate momentum regime the magnitude of the dif- 
ference is slightly less. It is somewhat remarkable that 
for both actions all three lattices agree in the ultraviolet. 
The tells us that ^(C, a) is essentially independent of a 
(but does differ between the FLIC Overlap and Wilson 
Overlap) . 



VI. CONCLUSION 

We have calculated the quark propagator in momen- 
tum space using the FLIC overlap fermion action in 
quenched QCD. A comparison with a previous study 
using the Wilson overlap Q indicates that both the 
mass function M(p) and the renormalisation function 
Z(q) agree well in the infrared regime. The compari- 
son was performed using a fitting technique to match the 
bare quark masses. The chiral behaviour of M(p) and 
Z(q) was studied, and clearly demonstrates non-linear 
behaviour at low momenta and light quark mass, which 
is consistent with expectations from chiral perturbation 



theory. 

The heavy quark mass behaviour of the FLIC overlap 
appears to be better than that of the Wilson overlap, as 
the FLIC overlap seems to inherit the benefits of smear- 
ing the irrelevant operators in the FLIC kernel. As one 
goes to lighter quark masses, M(p) for the two actions 
agree even in the ultraviolet. The same is not true of 
Z(q), as the two actions disagree in the ultraviolet for all 
masses. 

The scaling of the Wilson overlap and FLIC overlap 
results were also studied. The two actions showed sim- 
ilar scaling behaviour, with some differences. The re- 
sults for M(p) showed relatively good scaling across the 
three lattices at a = 0.096, 0.120 and 0.200. The scal- 
ing was best in the ultraviolet and in the very infrared, 
with some differences appearing in the intermediate mo- 
mentum regime. The results at a = 0.200 showed clear 
scaling violations at the heavy quark masses. Overall, 
for M(p) we see that the FLIC overlap shows superior 
scaling behaviour at heavy quark mass and in the ultra- 
violet. The scaling of the renormalisation function Z(p) 
was good in the ultraviolet, but disagreements appeared 
between the three lattices in the infrared. 

Given the superior properties of the FLIC overlap oper- 
ator outlined above, and the fact that the FLIC overlap 
operator is significantly cheaper to invert, we advocate 
the FLIC overlap operator as the operator of choice for 
chiral studies of hadron phenomenology. 

Future work will explore the effects of dynamical FLIC 
fermions j^] on the quark propagator [33]. Further in- 
vestigations of the potential effects of finite volumes and 
the chiral extrapolation function used will also be con- 
ducted. 
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FIG. 1: Cylinder cut data for the FLIC Overlap mass function M(p) at finite quark mass for the for the quenched lattices at 
(3 = 4.286 (top), P — 4.6 (middle) and /3 = 4.8 (bottom). The lefthand plots are against the discrete lattice momentum p, and 
the righthand plots are against the kinematical lattice momentum q. 
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FIG. 2: Cylinder cut data for the FLIC Overlap renormalization function Z^(p) (bottom) at finite quark mass for the quenched 
lattices at (3 = 4.286 (top), (3 — 4.6 (middle) and (3 = 4.8 (bottom). The lefthand plots are against the discrete lattice 
momentum p, and the righthand plots are against the kinematical lattice momentum q. The renormalization point for Zq (q) is 
£ = 10 GeV. The errant point in these plots is comes from the lightest quark mass. 
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FIG. 3: Cylinder cut data showing the dependence of the mass function M(p) for the FLIC Overlap (left) and the Wilson 
Overlap (right) on the bare mass m° at fixed momenta, for the lowest 10 momenta values. At small bare mass, the curves are 
ordered inversely to the momenta they represent, that is the smallest momenta is the topmost curve. Data is shown for for the 
quenched lattices at f3 = 4.286 (top), (3 = 4.6 (middle) and /3 = 4.8 (top). The solid curves are the quartic fits to the data. 
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FIG. 4: Cylinder cut data showing the dependence of the renormalisation function Z(p) for the FLIC Overlap (left) and the 
Wilson Overlap (right) on the bare mass m° at fixed momenta, for the lowest 10 momenta values. The curves are ordered 
according to the momenta they represent, that is the largest momenta is the topmost curve. Data is shown for for the quenched 
lattices at f3 = 4.286 (top), f3 = 4.6 (middle) and (3 — 4.8 (top). The solid curves are the quartic fits to the data. 
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FIG. 5: Cylinder cut data comparing the interpolated mass function M(p) (left) and renormalisation function Z(p) (right), for 
the three quenched lattices /3 = 4.286 (top), (5 = 4.60 (middle), /3 — 4.80 (bottom). The different curves from lowest to highest 
represent matched bare masses of m° = 0.0, 0.05, 0.1, 0.2, 0.4, 0.6 GeV. The solid points indicate the chiral limit. The different 
Z(q) have been offset vertically for clarity. 
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FIG. 6: Cylinder cut data showing the scaling of the FLIC overlap (top row) and Wilson overlap (bottom row) interpolated 
mass function M(p) (left) and renormalisation function Z(p) (right) for the three quenched lattices. The different curves from 
lowest to highest represent matched bare masses of m° = 0.0,0.1,0.2,0.4,0.6 GeV. The solid points indicate the chiral limit. 
The different Z(q) have been offset vertically for clarity. 



